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ABSTRACT 



Context. Close-in giant extrasolar planets ("Hot Jupiters") are believed to be strong emitters in the decametric radio 
range. 

Aims. We present the expected characteristics of the low-frequency magnetospheric radio emission of all currently known 
extrasolar planets, including the maximum emission frequency and the expected radio flux. We also discuss the escape 
of exoplanetary radio emission from the vicinity of its source, which imposes additional constraints on detectability. 
Methods. We compare the different predictions obtained with all four existing analytical models for all currently known 
exoplanets. We also take care to use realistic values for all input parameters. 

Results. The four different models for planetary radio emission lead to very different results. The largest fluxes are 
found for the magnetic energy model, followed by the CME model and the kinetic energy model (for which our results 
are found to be much less optimistic than those of previous studies). The unipolar interaction model does not predict 
any observable emission for the present exoplanet census. We also give estimates for the planetary magnetic dipole 
moment of all currently known extrasolar planets, which will be useful for other studies. 

Conclusions. Our results show that observations of exoplanetary radio emission are feasible, but that the number of 
promising targets is not very high. The catalog of targets will be particularly useful for current and future radio 
observation campaigns (e.g. with the VLA, GMRT, UTR-2 and with LOFAR). 
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1. Introduction 

In the solar system, all strongly magnetised planets are 
known to be intense nonthermal radio emitters. For a cer- 
tain class of extrasolar planets (the so-called Hot Jupiters), 
an analogous, but much more intense radio emission is ex- 
pected. In the recent past, such exoplanetary radio emission 
has become an active field of research, with both theoretical 
studies and ongoing observation campaigns. 

Recent theoretical studies have shown that a large va- 
riety of effects have to be considered, e.g. kinetic, magnetic 
and unipolar interaction between the star (or the stellar 
wind) and the planet, the influence of the stellar age, the 
potential role of stellar CMEs, and the influence of different 
stellar wind models. So far, there is no single publication 
in which all of these aspects are put together and where 
the different interaction models are compared extensively. 
We also discuss the escape of exoplanetary radio emission 
from its planetary system, which depends on the local stel- 
lar wind parameters. As will be shown, this is an additional 
constraint for detectability, making the emission from sev- 
eral planets impossible to observe. 

The first obser vation attempts go back at least to 
lYantis et alj (|1977l ). At the beginning, such observations 
were necessarily unguided ones, as exoplanets had not 



* Table 1 is only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http:/ /cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 



yet been discovered. Later observation campaigns con- 
centrated on known exoplanetary systems. So far, no 
detection has been achieved. A list and a comparison 
of past observati o n atte mpts can be found elsewhere 
(jGriefimeier et all . l2006al ). Concerning ongoing and fu- 
ture obser vations, studies a re perfor med or planned at 
the VLA (lLazio et all. |2004) . GMR T (iMaiid et all 120061: 
IWinterhalter et all I2006D. U TR2 (jRvabov et all |2004 ). 
and at LOFAR (|Farrell et all 12004 ) . To support these ob- 
servations and increase their efficiency, it is important to 
identify the most promising targets. 

The target selection for radio observations is based on 
theoretical estimates which aim at the prediction of the 
main characteristics of the exoplanetary radio emission. 
The two most important characteristics are the maximum 
frequency of the emission and the expected radio flux. The 
first p redictive studies (e.g. lZarka et allll997HFarrell et all 
1999) concentrated on only a few exoplanets. A first catalog 
containing estimations for rad io emission of a lar ge number 
of exoplanets was presented bv lLazio et al.l (|2004 ) . This cat- 
alog included 118 planets (i.e. those known as of 2003, July 
1) and considered radio emission energised by the kinetic 
energy of the stellar wind (i.e. the kinetic model, see below). 
Here, we present a much larger list of targets (i.e. 197 exo- 
planets found by radial velocity and/or transit se arches as 
of 2007, January 13, taken from http:/ /exopl anet.eu/[ ), and 
compare the results obtained by all four currently existing 
interaction models, not all of which were known at the time 
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of the previous overview study. As a byproduct of the ra- 
dio flux calculation, we obtain estimates for the planetary 
magnetic dipole moment of all currently known extrasolar 
planets. These values will be useful for other studies as, 
e.g., star-planet interaction or atmospheric shielding. 

To demonstrate which stellar and planetary parameters 
are required for the estimation of exoplanetary radio emis- 
sion, some theoretical results are briefly reviewed (section 
[2]) . Then, the sources for the different parameters (and their 
default values for the case where no measurements are avail- 
able) are presented (section [3]). In section |4] we present our 
estimations for exoplanetary radio emission. This section 
also includes estimates for planetary magnetic dipole mo- 
ments. Section [5] closes with a few concluding remarks. 



2. Exoplanetary radio emission theory 

2.1. Expected radio flux 

In principle, there are four different types of interaction be- 
tween a planetary obstacle and the ambient stellar wind, as 
both the stellar wind and the plane t can either be magne- 
tised or unmagnetised. [Zarkal (|2007l Table 1) show that for 
three of these four possible situations intense nonthermal 
radio emission is possible. Only in the case of an unmagne- 
tised stellar wind interacting with an unmagnetised body 
no intense radio emission is possible. 

In those cases where strong emission is possible, the ex- 
pected radio flux depends on the source of available energy. 
In the last years, four different energy sources were sug- 
gested: a) In the first model, the input power Pi npu t into 
the magnetosphere is assumed to be proportional to the to- 
tal kinetic energy flux of the solar wind protons impacting 
on the magnetopause (Dcsch & Kaiser, 1984; Zarka et al 
19971: IFarrell et all . Il999l: IZarka et all l2001bl: iFarrell et al 



2004 ; lLazio et al 



2004; IStevensl l2005t iGriefimeier et al 



2005. l2006bL 120071 ) b) Similarly, the input power P input into 



the magnetosphere can be assumed to be proportional to 
the magnetic energy flux or electromag netic Poynt i ng flux 
of the interp lanet ary magnetic field (I Zarka et all l2001bl : 
IFarrell et all 120041 : IZarkaL I2004L 1200(1 |2007j)- From the data 
obtained in the solar system, it is not possible to distinguish 
which of these models is more appropriate (the constants 
of proportionality i mplied in the relati ons given below are 
not well known, see lZarka et al.1 . [2001bT ). so that both mod- 
els have to be considered, c) For unmagnetised or weakly 
magnetised planets, one may apply the unipolar interaction 
model. In this model, the star-planet system can be seen 
as a g iant an alog to the Jupit er-Io system (| Zarka et all 
1200 It* IZarkaL I2004L I2006L 120071 ). Technically, this model is 
very similar to the magnetic energy model, but the source 
location is very different: Whereas in the kinetic and in the 
magnetic model, the emission is generated near the planet, 
in the unipolar interaction case a large-scale current sys- 
tem is generated and the radio emission is generated in the 
stellar wind between the star and the planet. Thus, the 
emission can originate from a location close to the stellar 
surface, close to the planetary surface, or at any point be- 
tween the two. This is possible in those cases where the solar 
wind speed is lower than the Alfve n veloc ity (i.e. for close- 
in planets, see e.g. iPreusse et all I2005T ) . Previous stud- 
ied have indicated that this emission is unlikely to be de- 
tectable, e xcept for sta r s with an ex t remely strong mag - 
netic field ijZarka et all l2001bL 12004 IZarkaL l200fit l2007h . 



Nevertheless, we will check whether this type of emission is 
possible for the known exoplanets. d) The fourth possible 
energy source is based on the fact that close-in exoplanets 
are expected to be subject to frequent and violent stel- 
lar eruptions (Khodachcn ko et all l2007bl ) similar to solar 
coronal mass ejections (CMEs). As a variant to the kinetic 
energy model, the CME model assumes that the energy for 
the most intense planetary radio emission is provided by 
CMEs. During periods of such CME-driven radio activity, 
considerably higher radio flux le vels can be achieved than 
durin g quiet stellar conditions (jGrieBmeier et all l2006bl 
120071) . For this reason, this model is treated separately. 

For t he kinetic energy case, the input power was first de- 
rived bv lDesch fc Kaiser ((1984), who found that it is given 

b y „ . 



-Pinput, kin OC nV c gR s . 



(1) 



In eq. JT]), n is the stellar wind density at the planetary 
orbit, w e ff is the velocity of the stellar wind in the reference 
frame of the planet (i.e. including the aberration due to 
the orbital velocity of the planet, which is not negligible 
for close-in planets), and R s denotes the magnetospheric 
standoff distance. 

The magnetic energy case was first discussed by 
IZarka et al.l (|2001tJ) . Here, the input power is given by 



Pinput, mag OC V c ffB j^R s 



(2) 



In eq. v e g is the velocity of the stellar wind in the 
reference frame of the planet, Bj_ if the component of the 
interplanetary magnetic field (IMF) perpendicular to the 
stellar wind flow in the reference frame of the planet, and 
R s denotes the magnetospheric stand off distance. 

For the unipolar interaction case (jZarka et alll2001bD . 
the input power is given by 



Pinput, unipolar OC V c ffB> j^R u 



(3) 



Eq. ([3]) is identical to eq. ([2]) , except that the obstacle is not 
the planetary magnetosphere, but its ionosphere, so that R s 
is replaced by Pi on , the radius of the planetary ionosphere. 

CME-driven radio e mission was first calculated by 
GrieBme ier et al.l (|2006b[ ). In that case, the input power is 
given by 



Pinput, kin, CME OC ncME%,OME^s 



(4) 



Eq. ([4]) is identical to eq. ([I]), except that the stellar wind 
density and velocity are replaced by the corresponding val- 
ues encountered by the planet during a CME. 

A certain fraction e of the input power Pi npu t given by 
eq. (QJ, ©, © or (LH} is thought to be dissipated within 
the magnetosphere: 



Pd = eP 



input 



(5) 



Observational evidence suggests that the amount of 
power emitted by radio waves P ra d is roughly p roportional 
to the power input Pinput (see, e.g. IZarkall2007l Figure 6). 
This can be written as: 



Pradio — pradio P<i — pradio ^Pi 



input 



(6) 



As Pd cannot be measured directly, one correlates the ob- 
served values of P ra dio with the calculated (model depen- 
dent) values of Pi npu t- Thus, one replaces Pi npu t by P ra dio on 
the left-hand side of the proportionalities given by ([T]), $2$, 
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^ and ((J) . The proportionality constant is determined by 
comparison with Jupiter. The analysis of the jovian radio 
emission allows to define three values for the typical radio 
spectrum: (a) the power during average conditions, (b) the 
average pow er during pe r iods o f high activity, and (c) the 
peak power (jZarka et al.l . [20041 . In this work, we will use 
the average power during periods of high activity as a ref- 
erence value for all four cases, with P ra dio„j = 2.1 • 10 11 W. 

The radio flux $ seen by an observer at a distance s from 
the emitter is rela t ed to the emitted radio power P ra dio by 
(|Griefimeier et all l2007h : 



^radio 



it 



(7) 



Here, SI is the solid angle of the beam of th e emitted 
radiation (Jl = 1.6 sr, see IZarka et all |2004(), and Af 
is the bandwidth of th e emission. We use A/ = /™ ax 
dGriefi mcicr etall l200l . where / c max is the maximum cy- 
clotron frequency. Depending on the model, P ra dio is given 
by eq. |T]), @, © or (HJ). The maximum cyclotron fre- 
quency /" lax is determined by the maximum magnetic field 
strength B™ ax close to the polar cloud tops (jFarrell et all 
11999ft : 



fc 



eBl 



2irm e 



enoM 
47r 2 r 



Ft 3 



24 MHz 



M 

~3 



(8) 



Here, m e and e are the electron mass and charge, R p is the 
planetary radius, /io is the magnetic permeability of the 
vacuum, and A4 is the planetary magnetic dipole moment. 
M. and R p denote the planetary magnetic moment and its 

radius relative to the respective valu e for Jupiter e.g. y Vf = 
M/Mj, with Mj = 1.56- 10 27 Am 2 dCain et al.l . fl995h and 
Rj = 71492 km. 

The radio flux expected for the four different models 
according to eqs. ([T]) to {7J and the maximum emission 
frequency according to fl5|) are calculated in section 0] for 
all known exoplanets. 

2.2. Escape of radio emission 

To allow an observation of exoplanetary radio emission, it 
is not sufficient to have a high enough emission power at 
the source and emission in an observable frequency range. 
As an additional requirement, it has to be checked that the 
emission can propagate from the source to the observer. 
This is not the case if the emission is absorbed or trapped 
(e.g. in the stellar wind in the vicinity of the radio-source), 
which happens whenever the plasma frequency 



plasma 



1 



2ir y eoTO e 



(9) 



is higher than the emission frequency at any point between 
the source and the observer. Thus, the condition of observ- 
ability is 

(10) 



jmiax max 
J plasma J c ' 



where /™ ax is taken at the radio source (e.g. the planet), 
whereas f p ™sma ^ s evaluated along the line of sight. As the 
density of the electrons in the stellar wind n decreases with 
the distance to the star, this condition is more restrictive 



at the orbital distance of the planet than further out. Thus, 
it is sufficient to check whether condition (flTJ|) is satisfied 
at the location of the radio-source (i.e. for n = n(d), where 
d is the distance from the star to the radio-source). In that 
case, the emission can escape from the planetary system 
and reach distant observers. In section [4] condition (fTO]) is 
checked for all known exoplanets at their orbital distance. 

Note however that, depending on the line of sight, not 
all observers will be able to see the planetary emission at all 
times. For example, the observation of a secondary transit 
implies that the line of sight passes very close to the plane- 
tary host star, where the plasma density is much higher. For 
this reason, some parts of the orbit may be unobservable 
even for planets for which eq. (flU)) is satisfied. 



2.3. Radiation emission in the unipolar interaction model 

An additional constraint arises because certain condi- 
tions are necessary for the generation of radio emis- 
sion. Planetary radio emission is caused by the cyclotron 
maser instability (CMI). This mechanism is only efficient 
in regions where the ratio between the electron plasma 
frequency and the electron cyclotron frequency is small 
enough. This condition can be written as 



/plasma ^ r\ a 

~ir~ ' 



(ii) 



where the electron cyclotron frequency f c is defined by the 
local magnetic field 



fc = 



iB 



(12) 



and /plasma is given by eq. ©. Observations seem to favor 
a critical frequency ratio close to the 0.1, while theoret- 
i cal work supports a critical f r equency ratio close to 0. 4 
(jLe Queau et all 1 19851 iHilgersl . 119921 : IZarka et all l2001al ). 
Fundamental O mode or second harmonic O and X mode 
emission are possible a lso for larger frequency ratios , but 
are much less efficient (jTreumannl . l2000t IZarkal , 120071) . To 
avoid ruling out potential emission, we use the largest pos- 
sible frequency ratio, i.e. f&y™ < 0.4. 

The condition imposed by eq. (fTTj) has to be fulfilled 
for any of the four models presented in section 12.11 For 
the three models where the radio emission is generated di- 
rectly in the planetary magnetosphere, n decreases much 
faster with distance to the planetary surface than B, so 
that eq. pip can always be fulfilled. For the unipolar inter- 
action model, the emission takes place in the stellar wind, 
and the electron density n can be obtained from the model 
of the stellar wind. In this case, it is not a priori clear 
where radio emission is possible. It could be generated any- 
where between the star and the planet. In section^ we will 
check separately for each planetary system whether unipo- 
lar interaction satisfying eq. (fTTj) is possible at any location 
between the stellar surface and the planetary orbit. 



3. Required parameters 

In the previous section, it has been shown that the de- 
tectability of planetary radio emission depends on a few 
planetary parameters: 

— the planetary radius R p 
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— the planetary magnetic moment M 

— the size of the planetary magnetosphere R s 

— the size of the planetary ionosphere i?i on 

— the stellar wind density n and its velocity u e g 

— the stellar magnetic field (IMF) B± perpendicular to 
the stellar wind flow in the frame of the planet 

— the distance of the stellar system (to an earth-based 
observer) s 

— the solid angle of the beam of the emitted radiation 

The models used to infer the missing stellar and plan- 
etary quantities require the knowledge of a few additional 
planetary parameters. These are the following: 

— the planetary mass M p 

— the planetary radius R p 

— its orbital distance d 

— the planetary rotation rate u 

— the stellar magnetic field (IMF) components B r ,B§ 

— the stellar mass M* 

— the stellar radius i?+ 

— the stellar age t* 

In this section, we briefly describe how these each of 
these quantities can be obtained. 



3.1. Basic planetary parameters 

As a first step, basic planetary characteristics have to be 
evaluated: 



d, Worbit and s are directly taken from the Extrasolar 
Planets Encyclopaedia at http : / / exoplan et . eu| as 
well as the observed mass M b s and the orbital eccen- 
tricity e. Note that in most cases the observed mass is 
the "projected mass" of the planet, i.e. M b s = M p smi, 
where i is the angle of inclination of the planetary orbit 
with respect to the observer. 

For eccentric planets, we calculate the periastron from 
the semi-major axis d and the orbital eccentricity e: 
rfmin = d/(l — e). Thus, the results for planetary radio 
emission apply to the periastron. 

For most exoplanets, the planetary mass is not pre- 
cisely known. Instead, usually only the projected mass 
MpSini is accessible to measurements, where i is 
the inclination of the planetary orbit with respect 
to the observer. This projected mass is taken from 
http : //exoplanet . eu| and converted to the median 
value of the mass: median(M p ) = Af b s ■ median (-^-r) = 

M bs « 1.15 Mobs- 
For transiting planets, R p is taken from original publi- 
cations. For non-transiting planets, the planetary radius 
is R p not known. In this case, we estimate the planetary 
radius based on its mass M p and orbital distance d, as 
explained in appendix [XJ The radius of a "cold" planet 
of mass M p is given by 



R P (d 



(aMp 



1/3 



1 



2/3 



1.47i?j ■ 



'1/3 



1 



2/3 



(13) 

with a — 6.1 • 10 4 m 3 kg 1 (for a planet with the same 
composition as Jupiter) and M max = 3.16 Mj. Again, 

Mp and M max denote values relative to the respective 



value for Jupiter (using Mj = 1.9 • 10 27 kg) 
of an irradiated planet is then given by 



R P (d) 



R p (d) 



R p (d = oo) # p (d = oo) 



1 + 0.05 



The radius 



T 

1 cq 



( 14 ) 

where T eq is the equilibrium temperature of the plane- 
tary surface. The coefficients To and 7 depend on the 
planetary mass (see appendix lATj . 



3.2. Stellar wind model 

The stellar wind density n and velocity v c g encountered 
by a planet are key parameters defining the size of the 
magnetosphere and thus the energy flux available to create 
planetary radio emission. As these stellar wind parameters 
strongly depend on the stellar age, the expected radio flux is 
a functio n of the estimated age of t he ex oplanetary host star 
([Stevens! . 120051: IGriefimeier et all l2005h . At the same time 
it is known that at close distances the stellar wind velocity 
has not yet reached the value it has at larger orbital dis- 
tances. For this reason, a distance-dependent stellar wind 
models has to be used to avoid overestimating the expected 
plane tary radio flux (|Griefimeie E l2006t IGriefimeier et all 
12007ft . 

It was shown (Griefimeier et all 12007ft that for stel- 
lar ages > 0.7 Gyr, the radial dependence of the stellar 
wi nd prop e rties c an be described by the stellar wind model 
of iParkerl (11958), an d that the more complex model of 
IWeber fc David (|l967ft is not required. In the Parker model, 
the interplay between stellar gravitation and pressure gra- 
dients leads to a supersonic gas flow for sufficiently large 
substellar distances d. The free parameters are the coro- 
nal temperature and the stellar mass loss. They are in- 
directly chosen by setting the stellar wind conditions at 1 
AU. More details on the model can be found e l sewhe re (e.g. 
iMann et all fl999t IPreusseL l200l iGriefimeierl . 12006ft . 

The dependence of the stellar wind density n and veloc- 
ity and i> e ff on the age of the stellar system is based on ob- 
servations of astrospheric absorption features of stars with 
different ages. In the region between the astropause and the 
astrospheric bow shock (analogs to the heliopause and the 
hcliosphcric bow shock of the solar system), the partially 
ionized local interstellar medium (LISM) is heated and 
compressed. Through charge exchange processes, a popu- 
lation of neutral hydrogen atoms with high temperature is 
created. The characteristic Lya absorption (at 1216 A) of 
this population was detectable with the high-resolution ob- 
servations obtained by the Hubble Space Telescope (HST). 
The amount of absorption depends on the size of the as- 
trosphere, which is a function of the stellar wind character- 
istics. Comparing the measured absorption to that calcu- 
lated by hydrodynamic codes, these measurements allowed 
the first empirical estimation of the evolutio n of the stel- 
lar mass loss rate as a function of ste llar age (|Wood et all 
12001 IWool 12001 iWood et all [2005h . It should be noted, 
however, that the res ulting estimates are only valid for stel- 
lar ages > 0.7 Gvr (IWood et all . |2005| ). From these ob- 
servations, (|Wood et alll2005T) calculate the age-dependent 
density of the stellar wind under the assumption of an age- 
independent stellar wind velocity. This leads to strongly 
overe stimated stellar wind densities, especially for youn g 
stars (jGriefimeier et all 120051: iHolzwarth fc Jardinell2007ft . 
For this reason, we combine these results with the model for 
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the ag e-dependence of the stellar wind veloci ty of lNewkirkl 
(1980). One obtains (| Griefimeier et all . l2007t ): 



u(lAU,£) =v 1 + 



t 



-0.43 



The particle density can be determined to be 

, \ -1.86±0.6 

n(lAU,*) = n ( 1 + - 



(15) 



(16) 



with v = 3971 km/s, n a = 1.04 • 10 11 m~ 3 and r = 2.56 • 
10 7 yr- 

For planets at small orbital distances, the keplerian ve- 
locity of the planet moving around its star becomes com- 
parable to the radial stellar wind velocity. Thus, the in- 
teraction of the stellar wind with the planetary magneto- 
sphere should be calculated using the effective velocity of 
the stellar wind plasma relative to th e planet, wh i ch take s 
into account this "aberration effect" (jZarka et all 120 01b). 
For the small orbital distances relevant for Hot Jupiters, 
the p lanetary orbits are circul a r because of tidal d i ssipa- 
tion dGoldreich fe Soterl 119661 : iDobbs-Dixon et all . l2004t 
lHalbwachs et all 12005 ). For circular orbits, the orbital ve- 
locity Worbit is perpendicular to the stellar wind velocity v, 
and its value is given by Kepler's law. In the reference frame 
of the planet, the stellar wind velocity then is given by 



Vcft 



- . /„,2 



rbit 



+ V 1 . 



(17) 



Finally, for the magnetic energy case, the interplane- 
tary magnetic field (B r ,B<p) is required. At 1 AU, the av- 
erage field stre ngth of the interplanetary magnetic field i s 
B im{ w 3.5nT (jMariani fc Neubauerl . [l990t iProlssl [20041. 
According to the Parker stellar wind model (|Parkeii ll958). 
the radial component of the interplanetary magnetic field 
decreases as 



-Bimf,r(d) = B r 



(18) 



This was later confirmed by Hel ios measurements. One finds 
B r . pa 2.6 nT and d = 1 AU (jMariani fc Neubauerl Il990t 
IProlssl 120041) . At the same time, the azimuthal component 
Amf,v3 behaves as 



d 



(19) 



I2004h . The average value of B lm [,g vanishes (B im i.e ~ 0). 
From B un { r (d) and B- lra t^(d), the stellar magnetic field 
(IMF) B±(d) perpendicular to the ste llar wind flow in the 
frame of the planet can be calculated (|Zarkall2007[) : 



with 



and 



a — arctan 



B, 



iraf,<p 



f3 = arctan ^ ^orbit ^ 



(20) 
(21) 
(22) 



We obtain the stellar magnetic field B+ relative to 
the solar magnetic field Bq under the assumption that 
it is inversely proportional to the rotation per i od P * 
(|Collier Cameron fc Jiaiikel 119941 iGriefimeier et al.l . l2007f) : 



B G 



(23) 



where we use P Q = 25.5 d and take Bq = 1.435 • 10~ 4 T 
as t he reference magneti c field strength at the solar sur- 
face (H reusse et all l2005f) . The st ellar rotation p eriod P+ is 
calculated from the stellar age t (jNewkirki . Il98d( ): 



P* oc 1 + 



0.7 



(24) 



where the time constant r is g iven by r = 2.56 • 10 7 yr 
(calculated from iNewkirkl . Il980h . 

Note that first measurements of stellar magnetic fields 
for planet-hosting stars are just bec oming availab l e from 
the spectropolarimeter ESPaDOnS (jCatala et all I2007D . 
This will lead to an improved understanding of stellar mag- 
netic fields, making more accurate models possible in the 
future. 



3.3. Stellar CME model 

For the CME-driven radio emission, the stellar wind param- 
eters n and v e g are effectively replaced by the corresponding 
CME parameters tt-cme and t> e ff,CME, potentially leading to 
much more intense radio emission than t hose driven by the 
kinetic energy of t he steady ste llar wind ( Gricfim eierl 120061 
IGriefimeier etail I2006U l200l . 

These CME param eters are estimated by 
iKhodachenko et"ail (|2007bt ). who combine in-situ mea- 
surements near the sun (e.g. by Helios) with remote solar 
observation by SoHO. Two interpolated limiting cases are 
given, denoted as weak and strong CMEs, respectively. 
These two classes have a different dependence of the 
average density on the distance to the star d. In the 
following, these quantities will be labeled ^cme(^) an d 
n GME( d ), respectively. 

For weak CMEs, the density ng ME (d) behaves as 



"cmeM =ncME..o(d/d ) 



-2.3 



(25) 



where the density at do = 1 AU is given by tIome o = 

ng ME (d = d ) = 4.9 • 10 6 m~ 3 . 

For strong CMEs, IKhodachenko et all (|2007ri) find 



'CME 



(d) 



n CME,0 



(d/do) 



-3.0 



(26) 



with B vfi ~ 2.4 nT (jMariani fc Neubauerl . Il990fc IProlssl witn n 



AU. 



CME.O 



'CME 



(d = d ) = 7.1 ■ 10 6 m~ 3 , and d = 1 



As far as the CME velocity is concerned, one has to 
note that individual CMEs have very different velocities. 
However, the average CME velocity v is approximately in- 
dependent of the subsolar distance, and is similar for both 
types of CMEs: 



v cme = v cme = "cme « 500 km/ s. 



(27) 



Similarly to the steady stellar wind the CME velocity given 
by eq. (j2Tj) has to be corrected for the orbital motion of the 
planet: 



"eff,CME 



M*G 



+ v. 



CME- 



(28) 
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In addition to the density and the velocity, the tempera- 
ture of the plasma in a coronal mass ejection is required for 
th e calculation of the size of th e magnetosphere. According 
to lKhodachenko et alJ (|2007bl lah. the front region of a CME 
consists of hot, coronal material (T ps 2 MK). This region 
may either be followed by relatively cool prominence mate- 
rial (T ps 8000 K), or by hot flare material (T ps 10 MK). 
In the following, the temperature of the leading region of 
the CME will be used, i.e. T CM e = 2 MK. 



3.4. Planetary magnetic moment and magnetosphere 

For each planet, the value of the planetary magnetic mo- 
ment M. is estimated by taking the geometrical mean 
of the maximum and minimum result obtained by differ- 
en t scaling laws. The assoc iated uncertainty was discussed 
by I Griefimeier et alj (12007th The differ en t scaling laws are 
compa red, e.g., by Farrell et al.1 (|!999f P1. IGrieBmeier et al.1 
(|2004 and lGriefimeierl (|2006O n order to be able to apply 



these scaling laws, some assumptions on the planetary size 
and structure are required. The variables required in the 
scaling laws are r c (the radius of the dynamo region within 
the planet), p c (the density within this region), er (the con- 
ductivity within this region) and oj (the planetary rotation 
rate). 



The size of the planetary core r c and its density p c The den- 
sity profile within the planet p(r) is obtained by describ- 
ing the planet as a polytropic gas sphere, using th e solu- 
tion of the Lane-Emd en equation ( Chandrasckharl Il957t 
ISanchez-Laveeal I2004D : 



p(r) 



TrMp 

ml 



sin 



(29) 



The size of the planetary core r c is found by search- 
ing for the value of r where the density p(r) be- 
comes large en ough for the t r ansiti o n to the liquid- 
meta l lic phase dSanchez-Lavegal 12004 IGrieBmeier et all 
I2005t iGriefimeierl . l2006fh Thetransition was assumed to 
occur at a density of 700 kg/m 3 , which is consistent wit h 
the range of parameters given bv ISanchez-Laveeal (|2004) . 
For Jupiter, we obtain r c — 0.85 Rj. 

The average density in the dynamo region p c is then 
obtained by averaging the density p(r) over the range < 
r < r c . For Jupiter, we obtain r c ps 1800 kg m~ 3 . 



The planetary rotation rate u> Depending on the orbital dis- 
tance of the planet and the timescale for synchronous rota- 
tion T sync (which is derived in appendix iBj) . three cases can 
be distinguished: 

1. For planets at small enough distances for which the 
timescale for tidal locking is small (i.e. T sync < 100 
Myr), the rotation period is taken to be synchronised 
with the orbital period (u) = Wf ps a> rbit)) which is 
known from measurements. This case will be denoted by 
"TL" . Typically, this results in smaller rotation rates for 
tidally locked planets than for freely rotating planets. 



2. Planets with distances resulting in 100 Myr < r sync < 
10 Gyr may or may not be subject to tidal locking. This 
will, for example, depend on the exact age of the plan- 
etary system, which is typically in the order of a few 
Gyr. For this reason, we calculate the expected char- 
acteristics of these "potentially locked" planets twice: 
once with tidal locking (denoted by "(TL)") and once 
without tidal locking (denoted by "(FR)"). 

3. For planets far away from the central star, the timescale 
for tidal locking is very large. For planets with r sync > 
10 Gyr, the effect of tidal interaction can be neglected. 
In this case, the planetary rotation rate can be as- 
sumed to be equal to the initial rotation rate u>i, which 
is assumed to be equal to the current rotation rate of 
Jupiter, i.e. w = wj with cjj = 1.77- 10~ 4 s _1 . This case 
will be denoted by "FR" . 



Note that tidal interaction does not perfectly synchro- 
nise the planetary rotation to its orbit. Thermal atmo- 
spheric tides resulting from stellar hea ting can drive planets 

away from synchron ous rota tion (IShowman & Guillot, 

120021 : ICorre a et all. 12003 iLaskar fcCorreial . l2004h. 
According to IShowman fc Guillotl (120021) . the correspond- 
ing error for u> could be as large as a f actor of two. On 
the ba sis of the example of r Bootis b, GrieB meier et al.l 
(|2007l ) show that the effect of imperfect tidal locking 
(in combination with the spread of the results found by 
different scaling laws) can lead to magnetic moments and 
thus emission frequencies up to a factor 2.5 higher than 
for the nominal case. Keeping this "error bar" in mind, we 
will nevertheless consider only the reference case in this 
work. 



The conductivity in the planetary core er c Finally, the con- 
ductivity in the dynamo region of extrasolar p lanets re- 
mains to be evaluated. According to Ncllis (2000), the elec- 
trical conductivity remains constant throughout the metal- 
lic region. For this reason, it is not necessary to average over 
the volume of the conducting region. As the magnetic mo- 
ment scaling is applied relative to Jupiter, only the relative 
value of the conductivity, i.e. cr/aj is required. In this work, 
the conductivity is assumed to be the same for extrasolar 
gas giants as for Jupiter, i.e. <j/<tj = 1. 



The size of the magnetosphere The size of the planetary 
magnetosphere R s is calculated with the parameters deter- 
mined above for the stellar wind and the planetary mag- 
netic moment. For a given planetary orbital distance d, 
of the different pressure contributions only the magneto- 
spheric magnetic pressure is a function of the distance to 
the planet. Thus, the st andoff distance R s is fo und from 
the pressure equilibrium (jGrieBmeier et all l2005h : 



Rs(d) 



8tt 2 (mn(d)v cS {d) 2 +2n(d)k B T) 



1/6 



(30) 



1 containing a typo in the sixth equation of the appendix. 



Here, m is the proton's mass, and /o is the form factor of 
the magnetosphere. It describes the magnetic field created 
by the magnetopause currents. For a realistic m agnetopause 
shape, a factor f Q = 1.16 is used (Voigt, 1995). In term of 
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units normalized to Jupiter's units, this is equivalent to 

1/6 



R s {d) ps 40i?j 



M 2 



n(d)v e s(dy 



1 n(d)k B T 



(31) 



with n(d) = n{d)/n } , v eS (d) = v cS (d)/v cS j, nj = 2.0 • 10 5 
m -3 , and v e s,i = 520 km/s. 

Note that in a few cases, especially for planets with very 
weak magnetic moments and/or subject to dense and fast 
stellar winds of young stars, eq. (f3"TJ)) yields standoff dis- 
tances R s < R p , where i? p is the planetary radius. Because 
the magnetosphere cannot be compressed to sizes smaller 
than the planetary radius, we set R s 



R p in those cases. 



3.5. Additional parameters 

The other required parameters are obtained from the fol- 
lowing sources: 



— The stellar ages t ± are taken from ISaffe et all (120051 
who g ave ages for 112 exoplanet host stars). ISaffe et alj 
(2005) compare stellar age estimations based on five dif- 
ferent methods (and for one of them they use two dif- 
ferent calibrations), some of which give more reliable 
results than others. Because not all of these methods 
are applicable to all stars, we use the age estimations 
in the following order or preference: chromospheric age 
for ages below 5.6 Gyr (using the D93 calibration), 
isochrone age, chromospheric age for ages above 5.6 
Gyr (using the D93 calibration), metallicity age. Note 
that error estimates for the two most reliable methods, 
namely isochrone and chromospheric ages, are already 
relatively large (30%-50%). In those cases where the age 
is not known, a default value of 5. 2 Gyr is used (th e 
median chromospheric age found by ISaffe et aLl . [2005) . 
This relatively high average age is due to a selection ef- 
fect (planet detection by radial velocity method is easier 
to achieve for older, more slowly rotating stars). As the 
uncertainty of the radio fl ux estimation be c omes very 
large for low stellar ages (|Griefi meier et all l2005h . we 
use a minimum stellar age of 0.5 Gyr. 

— For the solid angle of the beam, we assume the emis- 
sion to be analogous to the dominating contribu- 
tions of Jup i ter's radio emission and use = 1.6 sr 
(jZarkaet al.l . l2004D . 

4. Expected radio flux for know exoplanets 

4.1. The list of known exoplanets 

Table 1 shows what radio emission we expect from the 
presently known exoplanets (13.1.2007). It contains the 
maximum emission frequency according to eq. |5]l. the 
plasma frequency in the stellar wind at the planetary lo- 
cation according to eq. © and the expected radio flux ac- 
cording to the magnetic model {!]), the kinetic model 
and the kinetic CME model (U|). The unipolar interaction 
model is discussed in the text below. Table 1 also contains 
values for the expected planetary mass M p , its radius R p 
and its planetary magnetic dipole moment M. For each 
planet, we note whether tidal locking should be expected. 



Table 1 is only available in electronic form at the CDS. 



Note that |http : //exoplanet . eu contains a few more 
planets than table 1, because for some planets, essential 
data required for the radio flux estimation are not available 
(typically s, the distance to the observer). 

The numbers given in table 1 are not accurate results, 
but should be regarded as refined estimations intended to 
guide observations. Still, the errors and uncertainties in- 
volved in these estimation s can be considerable. As was 
shown in Griefi meier et al.l (|2007t ). the uncertainty on the 
radio flux at Earth, is dominated by the uncertainty 
in the stel lar age t+ ( for which the error is estimated as 
ps 50% by ISaffe et all l2005h . For the maximum emission 
frequency, /™ ax , the error is determined by the uncertainty 
in the planetary magnetic moment M, which is uncertain 
by a factor of two. For the planet r Bootes b, these effects 
translate into an uncertainty of almost one order of mag- 
nitude for the flux (the error is smaller for planets around 
stars of solar age), and an uncertainty of a factor of 2-3 
for the maximum emission frequency. T his error estimate 
i s deri ved and discussed in more detail bv lGriefimeier et al.l 
d2007h . 

The results given in table 1 cover the following range: 

— The maximum emission frequency found to lie beween 
to almost 200 MHz. However, all planets with / c max > 70 
MHz have negligible flux. Also note that any emission 
with / c < 5 to 10 MHz will not be detectable on Earth 
because it cannot propagate through the Earth's iono- 
sphere ("ionospheric cutoff"). For this reason, the most 
appropriate frequency window for radio observations 
seems to be between 10 and 70 MHz. 

— The radio flux according to the magnetic energy model, 
$sw,mag lies between and 5 Jy (for GJ 436 b). For 15 
candidates, $ S w,mag is larger than 100 mJy, and for 37 
candidates it is above 10 mJy. 

— The flux prediction according to the kinetic energy 
model, $sw,kin is much lower than the flux according to 
the magnetic energy model. Only in one case it exceeds 
the value of 10 mJy. 

— The increased stellar wind density and velocity during 
a CME leads to a strong increase of the radio flux when 
compared to the kinetic energy model. Correspondingly, 
^CME.kin exceeds 100 mJy in 3 cases and 10 mJy in 11 
cases. 

— Table 1 does not contain flux estimations for the unipo- 
lar interaction model. The reason is that the condition 
given by eq. ifTTj) is not satisfied in any o f the s t udied - 
cases . This is consistent with the result of IZarkal (2006, 
l2007h . who found that stars 100 times as strongly mag- 
netised as the sun are required for this type of emission. 
The approach taken in section 13.21 for the estimation 
stellar magnetic fields does not yield such strong mag- 
netic fields for stars with ages > 0.5 Gyr. Stronger mag- 
netic fields are possible (e.g. for younger stars). Strongly 
magnetised stars (even those without known planets) 
could be defined as specific targets to test this model. 

— The plasma frequency in the stellar wind, / PiSW , is neg- 
ligibly small in most cases. For a few planets, however, 
it is of the same order of magnitude as the maximum 
emission frequency. In these cases, the condition given 
by eq. (|10p makes the escape of the radio emission from 
its source towards the observer impossible. Of the 197 
planets of the current census, eq. (p~0|) is violated in 8 
cases. If one takes into account the uncertainty of the 
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stellar age (30-50%, see iSaffe et all f2005h . an uncer- 
tainty of similar size is introduced for the plasma fre- 
quency: In the example of r Bootis, the age uncertainty 
translates into a variation of up to 50% for the plasma 
frequency. With such error bars, between 6 and 14 plan- 
ets are affected by eq. (fTU|) . None of the best targets are 
affected. 

— The expected planetary magnetic dipole moments lie 
between and 5.5 times the magnetic moment of 
Jupiter. However, the highest values are found only for 
very massive planets: Planets with masses M < 2Mj 
have magnetic moments M < 2A4j. For planets with 
masses of M < Mj, the models predict magnetic mo- 
ments M < l.lMj. 

The results of table 1 confirm that the different mod- 
els for planetary radio emission lead to very different re- 
sults. The largest fluxes are found for the magnetic energy 
model, followed by the CME model and the kinetic en- 
ergy model. This is consistent with previous expectations 
(IZarka et all l2001bl: IZarkal 120061: iGrieflmeier etall l2006bt 
IZarkall200'7t iGriefimeier et all 120071 ). The univolar inter- 
action model does not lead to observable emission for the 
presently known exoplanets. Furthermore, the impact of 
tidal locking is clearly visible in the results. As it is cur- 
rently not clear which of these models best describes the 
auroral radio emission, it is not sufficient to restrict oneself 
to one scaling law (e.g. the one yielding the largest radio 
flux). For this reason, all possible models have to be consid- 
ered. Once exoplanetary radio emission is detected, obser- 
vations will be used to constrain and improve the models. 

Table 1 also shows that planets subject to tidal locking 
have a smaller magnetic moment and thus a lower maxi- 
mum emission frequency than freely rotating planets. The 
reduced bandwidth of the emission can lead to an increase 
of the radio flux, but frequently emission is limited to fre- 
quencies not observable on earth (i.e. below the ionospheric 
cutoff). 

The results of table 1 are visualized in figure Q] (for the 
magnetic energy model), figure^] (for the CME model) and 
figure [3] (for the kinetic energy model). The predicted plan- 
etary radio emission is denoted by open triangles (two for 
each "potentially locked" planet, otherwise one per planet). 
The typical uncertainties (approx. one order of magnitude 
for the flux, and a factor of 2-3 for the maximum emis- 
sion frequency) are indicated by the arrows in the upper 
right corner. The sensitivity limit of previous observation 
attempts are shown as filled symbols and as solid lines 
(a more detailed comparison of these observations can be 
found in IZarkal . [200l IGriefimeier et all 1200^ iGriefimeieR 
l2006f) . The expected sensitivity of new and future detec- 
tors (for 1 hour integration and 4 MHz bandwidth, or any 
equivalent combination) is shown for comparison. Dashed 
line: upgraded UTR-2, dash-dotted lines: low band and high 
band of LOFAR, left dotted line: LWA, right dotted line: 
SKA. The instruments' sensitivities are defined by the radio 
sky background. For a given instrument, a planet is observ- 
able if it is located either above the instrument's symbol or 
above and to its right. Again, large differences in expected 
flux densities are apparent between the different models. 
On average, the magnetic energy model yields the largest 
flux densities, and the kinetic energy model yields the low- 
est values. Depending on the model, between one and three 
planets are likely to be observable using the upgraded sys- 



tem of UTR-2. Somewhat higher numbers are found for 
LOFAR. Considering the uncertainties mentioned above, 
these numbers should not be taken literally, but should 
be seen as an indicator that while observation seem fea- 
sible, the number of suitable candidates is rather low. It 
can be seen that the maximum emission frequency of many 
planets lies below the ionospheric cutoff frequency, mak- 
ing earth-based observation of these planets impossible. A 
moon-based radio telescope however would give access t o 
radio emission with frequencies of a few MHz (Zarka, 2007). 
As can be seen in figures [TJ [3] and [3J this frequency range in- 
cludes a significant number of potential target planets with 
relatively high flux densities. 

Figures [1] [2] and [3] also show that the relatively high fre- 
quencies of the LOFAR high band and of the SKA telescope 
are probably not very well suited for the search for exoplan- 
etary radio emission. These instruments could, however, be 
used to search for radio emission generated by unipolar in- 
teraction between planets and strongly magnetised stars. 




1 MHz 10 MHz 100 MHz 1 GHz 

/ 

Fig. 1. Maximum emission frequency and expected radio 
flux for known extrasolar planets according to the magnetic 
energy model, compared to the limits of past and planned 
observation attempts. Open triangles: Predictions for plan- 
ets. Solid lines and filled circles: Previous observation at- 
tempts at the UTR-2 (solid lines), at Clark Lake (filled tri- 
angle), at the VLA (filled circles), and at the GMRT (filled 
rectangle) . For comparison, the expected sensitivity of new 
detectors is shown: upgraded UTR-2 (dashed line) , LOFAR 
(dash-dotted lines, one for the low band and one for the high 
band antenna) , LWA (left dotted line) and SKA (right dot- 
ted line). Frequencies below 10 MHz are not observable 
from the ground (ionospheric cutoff). Typical uncertainties 
are indicated by the arrows in the upper right corner. 



4.2. A few selected cases 

According to our analysis, the best candidates are: 

— HD 41004 B b, which is the best case in the magnetic 
energy model with emission above 1 MHz. Note that 
the mass of this object is higher than the upper limit 
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1 MHz 10 MHz 100 MHz 1 GHz 

/ 

Fig. 2. Maximum emission frequency and expected radio 
flux for known extrasolar planets according to the CME 
model, compared to the limits of past and planned obser- 
vation attempts. Open triangles: Predictions for planets. 
All other lines and symbols are as defined in figure [TJ 




1 GHz 



/ 

Fig. 3. Maximum emission frequency and expected radio 
flux for known extrasolar planets according to the kinetic 
energy model, compared to the limits of past and planned 
observation attempts. Open triangles: Predictions for plan- 
ets. All other lines and symbols are as defined in figure [TJ 



for planets (« 13A/j), so that it probably is a brown 
dwarf and not a planet. 

— Epsilon Eridani b, which is the best case in the kinetic 
energy model. 

— Tau Boo b, which is the best case in the magnetic energy 
model with emission above the ionospheric cutoff (10 
MHz). 

— HD 189733 b, which is the best case in both the mag- 
netic energy model and in the CME model which has 
emission above 1 MHz. 

— Gliese 876 c, which is the best case in the CME model 
with emission above the ionospheric cutoff (10 MHz). 



— HD 73256 b, which has emission above 100 mJy in the 
magnetic energy model and which is the second best 
planet in the kinetic energy model. 

— GJ 3021 b, which is the third best planet in the kinetic 
energy model. 

To this list, one should add the planets around Ups And 
(b, c and d) and of HD 179949 b, whose parent stars ex- 
hibit an increase of the ch r omospheric emission of about 
1-2% dShkolnik et all . I2TM [200l 120051) . The observations 
indicate one maximum per planetary orbit, a "Hot Spot" 
in the stellar chromosphere which is in phas e with the plan- 
etary orbit. The lead ang l es ob served by IShkolnik et al.l 
(|2003[ ) and IShkolnik et all (|2005l ) were recently explained 
with an Alfvcn-wing model using realistic stellar wind pa- 
rameters o btained from the stellar wind mode l by Weber 
and Davis (jPreussd . 120061 : iPreusse et all . |2006[) . This indi- 
cates that a magnetised planet is not required to describe 
the present data. The presence of a planetary magnetic field 
could, however, be proven by the existence of planetary ra- 
dio emission. Although our model does not predict high 
radio fluxes from these planets (see table 1), the high chro- 
mospheric flux shows that a strong interaction is taking 
place. As a possible solution of thi s prob l em, an inte nse 
stellar magnetic field was sugg ested (IZarkaL I2006ll2007h . In 
that case, table 1 underestimates the radio emission of Ups 
And b, c, d and HD 179949 b, making these planets inter- 
esting candidates for radio observations (eg. through the 
magnetic energy model or the unipolar interaction model) . 
For this reason, it would be especially interesting to obtain 
measurements of the stellar magnetic field for these two 
plane t-hosting stars (e.g. by the method of ICatala et all 
12001 . 

Considering the uncertainties mentioned above, it is im- 
portant not to limit observations attempts to these best 
cases. The estimated radio characteristics should only be 
used as a guide (e.g. for the target selection, or for statisti- 
cal analysis), but individual results should not be regarded 
as precise values. 

4.3. Statistical discussion 

It may seem surprising that so few good candidates are 
found among the 197 examined exoplanets. However, when 
one checks the li s t of cr iteria for "good" candidates (e.g. 
I Griefimeier et all l2006af ) , it is easily seen that only a few 
good targets can be expected: (a) the planet should be 
close to the Earth (otherwise the received flux is too weak). 
About 70% of the known exoplanets are located within 50 
pc, so that this is not a strong restriction, (b) A strongly 
magnetised system is required (especially for frequencies 
above the ionospheric cutoff). For this reason, the planet 
should be massive (as seen above, we find magnetic mo- 
ments M > 2A / lj only for planets with masses M > 2Aij). 
About 60% of the known exoplanets are at least as mas- 
sive as Jupiter (but only 40% have M v > 2.0Mj). (c) The 
planet should be located close to its host star to allow for 
strong interaction (dense stellar wind, strong stellar mag- 
netic field). Only 25% of the known exoplanets are located 
within 0.1 AU of their host star. 

By multiplying these probabilities, one finds that close 
(s < 50 pc), heavy (M p > 2.0Mj), close-in (d < 0.1 AU) 
planets would represent 8% (« 15 of 197 planets) of the cur- 
rent total if the probabilities for the three conditions were 
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independent. However, this is not the case. In the current 
census of exoplanets, a correlation between planetary mass 
and orbital distance is clea rly evident, with a lack of close- 
in massive planets (see e.g. lUdrv et"atl . I20Q3D . This is not a 
selection effect, as massive close-in planets should be easier 
to detect than low-mass planets. This correlation was ex- 
plained by the stronger tidal interaction effects for massive 
planets, leading to a faster decrease of the planetary orbital 
radius until the plane t reaches the stella r Roche limit and is 
effect ively destroyed (jPatzold &: Raueil l2002t Uiang et all 
|2003|) . Because of this mass-orbit correlation, the fraction 
of good candidates is somewhat lower (approx. 2%, namely 
3 of 197 planets: HD 41004 B b, Tau Boo b, and HD 162020 
b). 

4.4. Comparison to previous results 

A first comparative study of expected exoplanetary radio 
emission from a la rge number of planets was performed by 
lLazio et al.l ([2004), who compared expected radio fluxes of 
118 planets (i.e. those known as of 2003, July 1). Their 
results differ considerably from those given in table 1: 

— As far as the maximum emission frequency /™ ax is 
concerned, our results are consider a bly l ower than 
the frequencies given by lLazio et alj |2004). For Tau 
Bootes, their maximum emission frequency is six times 
larger than our result. For planets heavier than Tau 
Bootes, the discrepancy is even larger, reaching more 
than one order of magnitude for the very heavy cases 
(e.g. HD 168433c, for which they predict radio emis- 
sion with frequencies up to 2670 MHz). These differ - 
ence s have sever a l reas ons: Firstly, iFarrell et al.1 (|1999T ) 
and lLazio et all (|2004h assume that R p = Rj. Also, 
these wo r ks rel y on the magnetic moment scaling law of 
iBlackettl (jl947t ) , which has a large exponent in r c . This 
scaling law should not be u sed, as it was experimentally 
disproven (|Blackettl . 119521 ). Thirdly, these works make 

1 /3 

uses r c cx M p ' . Especially for planets with large masses 
like r Bootes, this yields unrealistically large core radii 
(even r c > R P in some cases), magnetic moments, and 
emission frequencies. Note that a good estimation of the 
emission frequency is particularly important, because a 
difference of a factor of a few can make the difference 
between radiosignals above and below the Earth's iono- 
spheric cutoff frequency. 

— The anticipated radio flux obtained with the kinetic en- 
ergy model < & g wjdn is much lower than the estimates of 
lLazio et al.1 (|2004f) . Typically, the difference is approxi- 
mately two orders of magnitude, but this varies strongly 
from case to case. For example, for Tau Boo b, our re- 
sult is smaller by a factor of 30 (where the difference is 
partially compensated by the low stellar age which in- 
creases our estimation), for Ups And b, the results differ 
by a factor of 220, and for Gliese 876 c, the difference is 
as large as a factor 6300 (which is partially due to the 
fact we take into account the small stellar radius and 
the high stellar age). 

— As was mentioned above, the analysis of the jovian radio 
emission allows to define three terms for the typical ra- 
dio spectrum: (a) the power during average conditions, 
(b) the average power during periods of high activity, 
and (c) the peak power (jZarkaet all 12004 . When com- 
paring the results of our table 1 to those of lLazio et al.l 



(|200l . one has to note that their table I gives the 
peak power, while the results in our table 1 were ob- 
tained using the aver age power du r ing p er iods of high 
activit y. Similarly to IFarrell et alj (|1999T ). lLazio et ail 
(|2004h assume that the peak power caused by varia- 
tions of the stellar wind velocity is two orders of magni- 
tud e higher than the a verage power. However, the val- 
ues |^rreir^t^aT| (|1999T ) use for average conditions corre- 
spond to periods of high activity, which are less than one 
order of magnitude below the peak power (jZarka et all 
2004) . During periods of peak emission, the value given 
in our table 1 would be increas ed by the same am ount 
(approximately a factor of 5, see lZarka et all [20041. For 
this reason, the peak radio flux is considerably overes- 
timated in these studies. 
— Estimated radio fluxes according to the magnetic energy 
model and the CME model have not yet been published 
for large numbers of planets. This is the first time the 
results from these models are compared for a large num- 
ber of planets. 

5. Conclusions 

Predictions concerning the radio emission from all presently 
known extrasolar planets were presented. The main param- 
eters related to such an emission were analyzed, namely 
the planetary magnetic dipole moments, the maximum fre- 
quency of the radio emission, the radio flux densities, and 
the possible escape of the radiation towards a remote ob- 
server. 

We compared the results obtained with various theo- 
retical models. Our results confirm that the four different 
models for planetary radio emission lead to very different 
results. As expected, the largest fluxes are found for the 
magnetic energy model, followed by the CME model and 
the kinetic energy model. The results obtained by the lat- 
ter model are found to be less optimistic than by previous 
studies. The unipolar interaction model does not lead to 
observable emission for any of the currently known plan- 
ets. As it is currently not clear which of these models best 
describes the auroral radio emission, it is not sufficient to 
restrict oneself to one scaling law (e.g. the one yielding the 
largest radio flux) . Once exoplanetary radio emission is de- 
tected, observations will be used to constrain and improve 
the model. 

These results will be particularly useful for the target 
selection of current and future radio observation campaigns 
(e.g. with the VLA, GMRT, UTR-2 and with LOFAR). We 
have shown that observation seem feasible, but that the 
number of suitable candidates is relatively low. The best 
candidates appear to be HD 41004 B b, Epsilon Eridani b, 
Tau Boo b, HD 189733 b, Gliese 876 c, HD 73256 b, and 
GJ 3021 b. The observation of some of these candidates is 
in progress. 
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Appendix A: An empirical mass-radius relation 

For the selection of targets for the search for radio emis- 
sion from extrasolar planets, an estimation of the expected 
radio flux <!> and of the maximum emission frequency /™ ax 
is required. For the calculation of these values, both the 
plane tary mass and the p l anetary radius are required (see , 
e.g. iFarrell et all Il999t iGrieBmeier et all . l2007t IZarkal 
2007). However, only for a few planets (i.e. the 16 presently 
known transiting planets) both mass and radius are known. 

In the absence of observational data, it is in principle 
possible to obtain planetary radii from numeric a l sim ula- 
tion, e.g. sim i lar to those of iBodenheimer et alj (|2003h or 
iBaraffe et all (I2003L 120051 ) . requiring one numerical run per 
planet. We chose instead to derive a simplified analytical 
fit to such numerical results. 



A.l. The accuracy of the fit 

The description presented here is necessarily only prelim- 
inary, as (a) numerical models are steadily further devel- 
oped and improved, and as (b) more transit observations 
(e.g. by the COROT satellite, which was launched recently) 
will provide a much better database in the future. This will 
considerably improve our understanding of the dependence 
of the planetary radius on various parameters as, e.g. plan- 
etary ma ss, orbital distance, or stellar metallicity (as sug- 
gested by iGuillot et all . l2006f ) . 

A. 1.1. What accuracy can we accept? 

Within the frame of the models presented in section [2 an 
increase in R p by 40% increases the expected radio flux 
by a factor of 2, and the estimated maximum emission fre- 
quency decreases by 40%. More generally, for a fixed plan- 
etary mass, 4> is roughly proportional to FlJ 3, and /™ ax is 
approximately proportional to Rp 1 . Thus, it appears that 
the assumption of a single standard radius for all plan- 
ets leads to a relatively large error. Comparing this to the 
other uncertainties involved in the estimation of rad i o char - 
acteristics (these are discussed in iGriefimeier et al"ll2007l) . 
it seems sufficient to estimate R p with 20% accuracy. 

A. 1.2. What accuracy can we expect? 

Several effects limit the precision in planetary radius we 
can hope to achieve: 

— The definition of the radius: The "transit radius" mea- 
sured for transiting planets is not exactly identical to 
the standard 1 bar radius. The differences are of the 
order of about 5-10%, but depend on the mass of the 



, but d ep 

planet (|Burrows et all . 120031 l2004f) . Because we com- 
pare modelled radii and observed radii without correct- 
ing for this effect, this limits the maximum precision we 
can potentially obtain. 

The abundance of heavy elements: The transiting planet 
around HP 149026 is sub stantially enriched in heavy el- 
ements (jSato et all 120051) . Models (|Bodenheimer et all 
|2003[ ) yield a smaller radius for planets with a heavy 
core than for coreless planets of the same mass (more 



than 10% difference for small planets). For a planet with 
unknown radius, a strong enrichment in heavy elements 
cannot be ruled out, as this case cannot be distinguished 
observationally from a pure hydrogen giant (i.e. one 
without heavy elements). 

For these reasons, we conclude that an analytical de- 
scription which agrees with the (numerical) data within 
~20% seems sufficient. For such a description, the error 
introduced by the fit will not be the dominant one. To get 
a better result, it is not sufficient to improve the approxi- 
mation for the radius estimation, but the more fundamental 
problems mentioned above have to be addressed. 

A. 2. An analytical mass-radius relation 

A. 2.1. The influence of the planetary mass 

A simple mass-radius relation, valid within a vast mass 
range , has been proposed by Lvnden-Bell k, O'Dwverl 
(|200lD and lLvnden-Bell fc Tout! (|2001l) : 



Rr, 



M, 



1/3 



(|ttpo 



1/3 



1 



2/3 



(A.l) 



The density po is depends on the planetary atomic compo- 
sition. Eq. (jA.lj) has a maximum in R p when M p — M max 
(corresponding to the planet of maximum radius). Fitting 
Jupiter, Sat urn and the plan et of maximum radius (R max — 
1.16.Rj, see iHubbardl . 119841) . we obtain M max = 3.16Afj 
and po — 394 kg m~ 3 for a Jupiter-like mixture of hydro- 
gen and helium (75% and 25% by mass, respectively). For 
a pure hydrogen planet, po = 345 kg m -3 . With there pa- 
rameters, eq. (|A.1I) fits w ell the results for cold planets of 
IBodenheimer et all (|2003l ). In this work, this estimation is 
used for the radii of non-irradiated (cold) planets. 

A. 2. 2. The influence of the planetary age 

It is known that the radius of a planet with a given mass 
depends on its age. According to m odels for the radii of 
isolated planets dBaraffe et all 120031 ). the assumption of a 
time-independent planetary radius leads to an error <11% 
for planets with ages above 0.5 Gyr and with masses above 
0.5Mj. In view of the uncertainties discussed above, this 
error seems acceptable for a first approximation, and we 
use 

Rp{M p ,t) « JZp(Afp). (A.2) 

A. 2. 3. The influence of the planetary orbital distance 

It is commonly expected that planets subjected to strong 
stellar radiation have a larger planetary radius than iso- 
lated, but otherwise identical planets. This situation is 
typical for "Hot Jupiters" , where strong stellar irradi- 
ation is suppo sed to dela y the planetary contraction 
(|Burrows et all 120001 120031 120041) . Clearly, this effect de- 
pends on the planetary distance to its star, d, and on the 
stellar luminosity L±. In the following, we denote the radius 
increase by r, which we define as 



R P (Mp,d) 



Rp(Mp,d- 



oo) 



(A.3) 
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Herein, R p (M p ,t, d) denotes the radius of a planet under 
the irradiation by its host star, and R p (M p , t,d — oo) is the 
radius of a non-irraditated, but otherwise identical planet. 

Different exoplanets have vastly different host stars. A 
difference of a factor of two in stellar mass can result in a 
difference of more than order of magnitude in stellar lumi- 
nosity. For this reason, it is not sufficient to take the orbital 
distance as the only parameter determining the radius in- 
crease by irradiation. 

Here, we select the equilibrium temperature of the 
lanetary surface as the basic parameter. It is defined as 
Bodenheimer et all . [20031 ) 



(1-A)L* 



167R7s B d 2 (l + e 2 /2) 2 



1/4 



(A.4) 



where the planetary albedo is set as A — 0.4 in the follow- 
ing. The stellar luminosity L+ is calculated fro m the stellar 
mass according to the analytical fit given by I Tout et all 
(1996) for the zero-age main sequence. In this, we as- 
sumed solar metallicity for all stars. Finally, <tsb denotes 

the Stefan-Boltzman n constant. 

Using the data of lBodenheimer et al.l (|2003h . we use the 
following fit for the radius increase r: 



' rp \ 7 

(..I).', 

Jo 



(A.5) 



This form was selected because it has the correct qualita- 
tive behaviour: It yields a monotonous decrease of r with 
decreasing equilibrium temperature T cq (increasing orbital 
distance d). In the limit T cq — > (d — > oo) we fi nd r — > 1 . As 
an alt ernative to a continuous fit like eq. (|A.5p . lLazio et all 
(2004) use a step-function, i.e. assume an increased ra- 
dius of r = 1.25 for planets closer than a certain distance 
d± = 0.1 AU only. However, this treatment does not re- 
duce the number of fit-constants. Also, for planets with 
orbital distances close to d\, the results obtained with a 
step-function strongly depend on the somehow arbitrary 
choice of d%. A continuous transition from "irradiated" to 

"isolated" planets is less prone to this effect. 

The numerical results of iBodenheimer et ail (|2003l ) 
show that the ratio r also depends on the mass of the planet: 
for small planets, r = 1.4 for an equilibrium temperature of 
2000 K, whereas for large planets, r < 1.10. For this reason, 
we allow to coefficients T and 7 to vary with M p : 



and 



7 = 1.15 + 0.05- 



(A.6) 



(A.7) 



With the set of coefficient c M = 764 K , c t , 2 = 0.28, 
C7 1 = 0.59Afj and c 7 2 = 1- 03, we obtain an analytica l 
fit to the numerical results of IBodenheimer et al.l (|2003h . 
The maximum deviation from the numerical results is be- 
low 10% (cf. Fig. EH)- 

For comparison, Fig. IA.ll also shows the value of r for 
the transiting exoplanets OGLE-TR-lOb, OGLE-TR-56b, 
OGLE-TR-lllb, OGLE-TR-113b, OGLE-TR-132b, XO- 
lb, HD 189733b, HD 209458b, TrES-lb and TrES-2b as 
small crosses. Here, r is calculated as the ratio of the ob- 
served value of R p (M p , d) and the value for R p (M p , d = 00) 



calculated according to eq. (jA.ip . As the mass of all tran- 
siting planets lies between O.llMj < M p < 3.0Mj, one 
should expect to find all crosses between the two limiting 
curves. For most planets, this is indeed the case: only for 
HD 209458b, r is considerably outside the area delimited 
by the two curves. Different explanations have been put 
forward for the anomalously large radius of this planet, but 
so far n o conclusive answer to this question has been found 
(see e.g. lGuillot et al.l . l2006l and references therein). As nu- 
merical models cannot reproduce the observed radius of this 
planet, one cannot expect our approach (which is based on 
a fit to numerical results) to reproduce it either. For the 
other planets, Fig. I A. II shows that our approach is a valid 
approximation. 



radius increase by irradiation 




Teff [K] 

Fig. A.l. Radius increase r as a function of the equilib- 
rium temperature T cq according to eq. (|A.5|) . (|A.6|1 and 
(IA.7[) . Upper line: radius increase r due to irradiation for 
planets of mass M p = O.llMj. Lower line: Upper line: 
radius increase r due to irradiation for planets of mass 
M p = 3.0 Mj. Open symbols: results of the numerical cal- 
culation of lBodenheimer et al.l (|2003h for planetary masses 
M p = O.llMj (diamonds) and M p = 3.0Mj (triangles). 
Crosses: radius increase for observed transiting planets (see 
text). 



Appendix B: Tidal locking? 

For the estimation of the planetary magnetic dipole mo- 
ment in section [3T4l we require the planetary rotation rate. 
This rotation rate rotation greatly depends on whether the 
planet can be considered as tidally locked, as freely rotat- 
ing, or as potentially locked. In this appendix, we discuss 
how we evaluate the tidal locking timescale r sync , which 
decides to which of these categories a planet belongs. 

The value of the planetary rotation u> depends on its dis- 
tance from the central star. For close-in planets, the plan- 
etary rotation rate is reduced by tidal dissipation. In this 
case, tidal interaction gradually slows down the planetary 
rotation from its initial value w; until it reaches the final 
value u>i after tidal locking is completed. 

It should be noted that for planets in eccentric orbits, 
tidal interaction does not lead to the synchronisation of 
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the planetary rotation with the orbital period. Instead, the 
rotation period also depends on the orbital eccentricity. At 
the same time, the t imescale to reach this eq uilibrium rota- 
tion rate is reduced (|Laskar fc Correial I2004D . Similarly, for 
planets in a n oblique orbit , the e quilibrium rotation period 
is modified (|Levrard et all [2007). However, the locking of a 
hot Jupiter in a non-synchronous spin-orbi t resonance ap- 
pears to be unlikely for distances < 0.1 AU (jLevrard et all 
2007). For this reason, we will calculate the timescale for 
tidal locking under the assumptions of circular orbits and 
zero obliquity. In the following, the tidal locking timescale 
for reaching uoi is calculated under the following simplifying 
assumptions: prograde orbit, spi n parallel to orbit (i.e. zero 
obliquity), and zero eccentricity ([Murray fc Dermottl [l999. 
Chapter 4). 

The rate of change of the planetary rotation velocity uj 
for a planet with a mass of M p and radius of R D around 
a star of mass is giv en by (|Goldreich fc Soteil Il966t 
iMurrav fc Dermottl . [l999h : 



dtu _ 9 1 (GM P \ (MA 2 fR p 



dm 

~dt ~ AaQ ' V Rl 



Mr) 



(B.l) 



where the constant a depends on the internal mass distri- 
bution within the planet. It is defined by a = I /(M p R p ), 
where I is the planetary moment of inertia. For a sphere of 
homogeneous density, a is equal to 2/5. For planets, gener- 
ally a < 2/5. Q' p is the modified Q- value of the planet. It 
can be expressed as (jMurrav fc Dermottl . Il999f ) 



Q' 



3Q P 



2k 



2,P 



(B.2) 



where fc 2 . p is the Love number of the planet. Q p is the plan- 
etary tidal dissipation factor (the larger it is, t he sm aller 
is the tidal dissipa t ion), defined by[MacDonald (1964j) and 
iGoldreich fc"Sotei1 (fl9Ml . 

The time scale for tidal locking is obtained by a com- 
parison of the planetary angular velocity and its rate of 
change: 

7"sync — . ■ v^ - "/ 

A planet with angular velocity u>\ at t = (i.e. after forma- 
tion) will gradually lose angular momentum, until the an- 



gular velocity reaches ojf at t = r sync . Insertion of eq. (jB.l 
into eq. (|B.3[) yields the following expression for r ( 



sync ■ 



Isync ~ g OL Q p 



R 3 P 



. (M p 



(B.4) 



The importance of this effect strongly depends on the 
distance (r sync cx d 6 ). Thus, a planet in a close-in orbit 
(d < 0.1 AU) around its central star is subject to strong 
tidal interaction, leading to gravitational locking on a very 
short timescale. 

In the following, we briefly describe the parameters (a, 
Q' p , uji and ojf) required to calculate the timescale for tidal 
locking of Hot Jupiters. 

Structure parameter a For large gaseous planets, the equa- 
tion of state can be approximated by a polytrope of index 
k = 1. In that case, the structure p arameter a (defi ned by 
a = I/MpRl) is given by a = 0.26 (|Gu et all l2003fl . 



Tidal dissipation factor Q' p For planets with masses of 
the order of one Jupiter mass, one finds that k2. P has a 
value of fc2. D ~ 0-5 for Jupiter (jMurrav fc Dermottl . 1999; 
i Laskar fc Correial 120041) and fe .p ~ 0.3 for Saturn ([Peak, 
Il999t lLaskar fc Correial j2004j).'The value of k 2 , p = 0.5 
will be used in this work. With eq. (|B.2[) . this results in 
Q P «3Q P . 

For Jupiter, one finds the f ollowing ran ge of allowed val- 
ues: 6.6 ■ 10 4 < Q p < 2.0 • 10 6 (lPeald . ll999D . Several estima- 
tions of the turbulent dissipation within Jupiter yield Up- 
values larger than this upper limit, while othe r theories pre- 
dict v alues consist ent with this upper limit (jMarcv et all 
119971 iPeald . 119991 and references therein). This demon- 
strates that the origin of the value of Q D is not well un- 
derstood even for Jupiter (jMarcv et all Il997l ) . 

Extrasolar giant planets are subject to strongly different 
conditions, and it is difficult to constrain Q p . Typically, Hot 
Jupiters are assumed to behave similarly to Jupiter, and 
values in the range 1.0- 10 5 < Q' p < 1.0- 10 6 are used. In the 
following, we will distinguish three different regimes: close 
planets (which are tidally locked), distant planets (which 
are freely rotating), and planets at intermediate distances 
(which are potentially tidally locked) . The borders between 
the "tidally locked" and the "potentially locked" regime is 
calculated by setting r sync = 100 Myr and Q' p = 10 6 . The 
border between the "potentially locked" and the "freely 
rotating" regime is calculated by setting r sync = 10 Gyr and 
Q' p = 10 5 . Thus, the area of "potentially locked" planets is 
increased. 



Initial rotation rate lo{ The initial rotation rate uii is not 
well constrained by planetary formation theories. The rela- 
tion between the planetary angular momentu m density and 
plane tary mass observed in the solar system (jMacDonaldl 
I1964D sug gests a pr i mordi al rotation period of the order of 
10 hours (jHubbardl . Il98l Chapter 4). We assume the ini- 
tial rotation rate to be equal to the current rotation rate of 
Jupiter (i.e. ui = u>i = ujj) with wj = 1.77 • 10~ 4 s . 



Final rotation rate uif As far as eq. (|B.4j) is concerned, Uf 
can be neglected (jGricfimcierl 120061 ). 
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